Sensitivity of the Jet Quenching Observables to 
the Temperature Dependence of the Energy Loss 
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The quenching of minijet (particles with pr » T, Aqcd ) in ultra-relativistic heavy-ion collisions 
has been one of the main prediction and discovery at RHIC. We analyze the correlation between 
different observables like the nuclear modification factor Raa(pt), the elliptic flow and the ratio of 
quark to gluon suppressions. We show that the temperature (or entropy density) dependence of the 
in- medium energy loss strongly affects the relation among these observables. In particular the large 
elliptic flow and the nearly equal Raa(pt) of quarks and gluons can be accounted for only if the 
energy loss occurs mainly around T c and the q o g conversion is significant. The use of an equation 
of state fitted to lattice QCD calculations, slowing down the cooling as T — > T c , seems to contribute 
to both the enhancement of V2 and the efficiency of the conversion mechanism. 

PACS numbers: 12.38.Mh , 25.75 -q , 25.75 Ld 



The experiments at the Relativistic Heavy Ion Col- 
lider (RHIC) are dedicated to study the properties of the 
matter at exceedingly high density and temperature. At 
such extreme conditions the matter was expected to un- 
dergo a deconfinement and chiral phase transition and 
have quark and gluons as degrees of freedom. The theo- 
retical and experimental efforts have shown that indeed 
a new form of matter has been created [1] . Such a matter 
appears as a nearly perfect fluid with very low viscosity 
to entropy density [2, 3]. It develops strong collective 
modes with quarks as degrees of freedom and hadronizcs 
in a modified way respect to pp collisions, at least in the 
intermediate 2 < px < 5 GeV region [4, 5]. One way to 
probe the created matter is to exploit the high energy 
jets (pr » T,Aqcd) produced by the hard collisions 
at the initial stage. They are internal probes propagat- 
ing through the fireball and interacting with the medium, 
hence carrying information on its properties as proposed 
long ago in Ref.s [6-8]. It has indeed been shown that 
the matter has a very high opacity respect to high px 
partons that traverse the hot medium in agreement with 
the expectations about the energy loss in QCD medium 
[8-10]. This energy loss can be quantified by the sup- 
pression of observed hadron spectra at high transverse 
momenta pt, as well as in the suppression of back-to- 
back di-hadron correlations with a high-p t trigger, when 
compared with pp or d-A collisions [1, 11]. Both these 
phenomena related to the "jet quenching" have been ob- 
served and represent one of the major discoveries of the 
RHIC experimental program [1]. 

However even if the observation of the jet suppression 
cannot be questionated there are several fundamental 
questions that still remain open. There are indeed several 
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models that depend in a different way on temperature T 
[12, 13], or that do not depend explicitly on T but on the 
q transport coefficient [14, 15], some others are based on 
perturbative approach other on higher twist expansion 
[16, 17]. 

Despite differences all the approaches seems to be able 
to describe the amount of suppression Raa{pt) observed 
experimentally. It is patent that one needs to go one or 
two steps further and fortunately experimentally there 
is already the availability of other observables related to 
the jet quenching phenomena. Interestingly the present 
models does not seem to be able to account for all of them 
simultaneously. In this paper we will focus on two observ- 
ables beyond the Raa(pt)- One is the elliptic flow ^(pt) 
and the other is the flavor dependence of the suppression 
that we will mainly discuss in terms of Raa{q) / 'Raa{q) ■ 
The latter can be experimentally inferred by a systematic 
comparison of the different suppression for n, p, K,p,p 
that are differently related to quark and gluon suppres- 
sion. 

The purpose of the present paper is to show that the 
jet-quenching mechanism carries much more information 
than what can be inferred from only the nuclear modifi- 
cation factor Raa(.Pt) and even the jet-triggered angular 
information. We suggest that the study of the correlation 
between the elliptic flow V2(pt) and the flavor depen- 
dence of the quenching Raa(<i) / ' Raa{9) is rich of infor- 
mation on the temperature dependence of the quenching 
and on the mechanism of parton flavor conversion. 

We point out that a correlation between t>2 (pt) and 
Raa(q) / RAA(g) is sensitive to the temperature (or en- 
tropy density s) dependence of AEi oss and on the den- 
sity profile of the bulk. The latter becomes unexpectedly 
dramatic if an extreme T, p or s dependence is consid- 
ered as in the recent works of E. Shuryak and J. Liao 
[18] or of V.S. Pantuev [19]. This shows in general that 
once one goes beyond the Raa (jPT ) a more careful treat- 
ment of the time evolution of the fireball becomes manda- 
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A[GeV] B[GeV] n f 
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1440 


1.5 


8.0 


q 


670 


1.6 


7.9 


q 


190 


1.9 


8.9 



TABLE I: Parameters for initial parton minijet distribu- 
tion given in Eq.(l) at midrapidity for Au+Au at ^snn — 
200AGeV. Taken from Ref. [20]. 



tory, but also the simultaneous description of Raa(pt), 
v 2 {pt), and Raa(q)/ Raa{q) ( or the ratio of the Raa(pt) 
between different hadrons) is non trivial and more rich 
of information. In particular we find that an energy loss 
that increase as T — > T c , the q o g in- medium conver- 
sion and the expansion-cooling of the fireball according 
to a lattice QCD EoS all go in the direction of improving 
the agreement with the three observables. 

The paper is organized in five Sections. In Section I the 
main ingredients of the model are presented. In Section 
II the model is applied to calculate the Raa(pt) as a 
function of momentum and centrality. In Section III the 
elliptic flow and the ratio of quark to gluon Raa(pt) are 
discussed together with their correlation. In Section IV 
the impact of a realistic equation of state is presented. In 
the last Section we summarize the conclusions from the 
present study. 



I. MODELLING THE JET QUENCHING 

Our modelling of the jet energy loss is based on the 
very widely used adiabatic-like approximation for which 
the jet loses energy in a bulk medium that is indepen- 
dently expanding and cooling. Therefore the jet energy 
loss is considered as a small perturbation of the bulk dy- 
namics. This is essentially the main assumption in the 
model till now developed, even if the higher energy at 
LHC should make this assumption inadeguate due to the 
significant amount of energy that is contained in the high 
Pt partons intially produced. 

The main ingredients of our model can be easily 
sketched in the following way. 

Initial Conditions - The parton distributions are calcu- 
lated in the NLO pQCD scheme: They are parametrized 
by power law functions as: 



dNf _ A f 

(1 +PT /Bf) n ' 



(1) 



with / = q,q,g. The transverse momentum px is in 
unit of GeV and the value of the parameters Af , B f ,rif 
is given in Table I and are taken from Ref. [20]. Such a 
choice is driven mainly by the intention to make a direct 
contact and comparison with the several works discussing 
jet quenching with the flavor conversion [20-24] . 



As regards the parton distribution in space coordi- 
nates, it scales with the number of binary nucleon col- 
lision N co ii according to the standard Glauber model 
[25, 26]. The initial conditions for the bulk medium are 
described by the density profile p(f, z, r) that in the lon- 
gitudinal direction evolves according to the Bjorken ex- 
pansion at the velocity of light. The initial transverse 
density profile is instead proportional to the standard 
Glauber model participant distribution: 



P P art(b,r) = t A {f) [l - e -— tA ^ 
+t B (b-r) \\- e -' 7NNtA ^ 



(2) 



with <tnn = 42mb, b the impact parameter vector and t A 
the nuclear thickness function normalized to the number 
of nucleons and given by: 
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(3) 



where p A (r) is the nuclear density that we have taken 
to be a Wood-Saxon (WS) with radius R = 6.38/m and 
thickness a = 0.535 /to. The total number of participants 
is therefore 



N par t{b) 



j d 2 rp part {b,r) (4) 



Beyond the Glauber density profile we have also con- 
sidered a simplified sharp elliptic (SE) shape with the x 
and y axis adjusted to reproduce the same eccentricity 
of the Glauber model at each impact parameter. The SE 
has been used for the description of the bulk in several 
jet quenching modeling [20-22, 27]. 

Bulk Density Evolution - For the sharp elliptic shape 
the density space-time evolution is given by: 



P{x,y,r) 



1 



dN 



tA t {t) dy k 



9 1 



R%{r) R?Jt) 



(5) 



where At = irRlR?, is the transverse area of the evolv- 
ing fireball and R x ,R y is the length of the two axes of 
the ellipse. These can in general evolve according to an 
expansion at constant acceleration or constant velocity: 



R x {t) = R xo + v t t + i(a T + e a )r 2 
Ry{r) = R yo + v t t + i(a T - e a )r 2 



(6) 



with vt the transverse velocity, ar±£ a the acceleration 
that through e a can be taken to be different between the 
x and y direction in order to simulate the anisotropic 
azimuthal expansion that reduces the eccentricity with 
time. We have used a typical value of a T — 0.08 fm -1 
that generates a final radial flow /3 = 0.4 with = 0, 
e a = 0.04 fm -1 with a slight dependence on the impact 
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parameter of the collision. However the sensitivity of the 
jet quenching on these parameters is quite limited. 

For the case with the Glauber profile the evolution of 
the local density is similarly given by 

where respect to Eq.(5) the thcta function is substituted 
by the profile function P e ff(x,y,r) and At is now the 
effective area given by the space integral of the profile 
function: 



P eff (x,y,T 



Npart(x,y,T) 



A T = 



part(0,0,r) 

J J dxdyP eff (x,y 7 T) 



(8) 



where N part {x, y, To) is given by the Glauber model while 
the time dependence is determined by the expansion in 
the (x,y) plane according to Eq.(6). 

Energy Loss - The aim of our work is to explore the 
consequences of different T dependences of energy loss on 
the observables, hence we employ various schemes for the 
energy loss. However to make a connection to the large 
amount of effort to evaluate gluon radiation in a pQCD 
frame we will use also the GLV formula at first order in 
the opacity expansion [10, 27]: 



A£(p,t, m ) 9tt 3 (2E 
— = —C R a s p{x,y,r)Tlog 



(9) 



where Cr is the Casimir factor equal to 4/3 for quarks 
and 3 for gluons, \i = gT is the screening mass with g = 3 
in agreement with 1QCD results [28] and r is the time mi- 
nus the initial time t% — 0.2 fm. Furthemore considering 
massless partons at midrapidity E — pt- There are cor- 
rections to Eq. (9) coming from higher order that can be 
approximately accounted for by a rescaling Z factor of the 
energy loss. However this not really relevant for the ob- 
jectives of the present work because we will renormalize 
the energy loss in order to have the observed amount of 
suppression Raa(pt) for central collisions, see next Sec- 
tion. In fact our purpose is to study how other variables 
can change once Raa(j>t) has been fixed to experimetal 
data for central collisions. 

Usually in the GLV, as well as in other approaches, 
the temperature evolution of the strong coupling a s is 
discarded. We will consider the impact of such a depen- 
dence to understand the amount of T dependence that 
can come simply from the asymptotic freedom in a pQCD 
approach. The scale dependence of the strong coupling 
can be writte as: 



a s (Q 2 



-in 



0o H-Qi/Aqcd) 



(10) 



where 0q = 11 — -Nf and the thermal scale Q 2 = 

(27rT) 2 which allows to get a correct behavior of the 
screening mass fi on the energy scale [28]. 
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FIG. 1: Temperature dependence of the energy loss for a 
parton at pr = 10 GeV. By dashed and dot-dashed line the 
GLV energy loss with constant and T dependent a a coupling 
(see text). The thick line is the case in which the energy loss 
take place only closer to the phase transition and the thin 
line represents an opposite case in which the energy loss take 
place only at high T. 



In Fig. 1 we show by dot-dashed and dashed lines the 
temperature dependence of the energy loss for the GLV 
with a dependence of the coupling according to Eq.(10) 
(GLV-cv s (T)) and with a constant coupling ct s = 0.27 
(GLVc). In itself the effect of the asymptotic freedom 
that reduces a s at increasing energy scale, significantly 
modifies the temperature dependence of the energy loss. 
However we will see that such effect is not very large once 
the coupling itself is readjusted to produce the correct 
amount of suppression. 

In Fig.l we have also shown two other opposite cases 
for AE/At: the thick line that shifts the energy loss to 
lower temperature (hence low density p or entropy den- 
sity s) as suggested in [18, 19] and the other (thin line) 
that gives a dominance of quenching at high T, consid- 
ered here just for comparison respect to the opposite case. 
For both Ei oss we have a standard dependence on the 
transverse momentum, AE(pt) ~ p T with 7 = 0.6. 

We have plotted in Fig.l the AEi oss as a function of 
the temperature, however under the hypothesis of local 
equilibrium for the bulk this is equivalent to a density or 
an entropy-density dependence used by other authors. In 
the model implemented we generally employ the free gas 
approximation to relate density, temperature and time. 
As well known one has: 



P 



C(3) /3 



T 



(11) 



where the last expression is obtained with d q> q = 24, d g = 
16. From Eq. (11) with an initial temperature To = 340 
MeV at To = 0.6 fm and a transverse area At ~ 90 fm 2 
one has a dN/dy ~ 1000 for b = 3/m corresponding to 
— 10% in agreement with standard estimates. Eq.(ll) 
allows one to evaluate the local temperature from the 
local density given by Eq.(5) for the SE profile and by 
Eq.(7) for the Glauber WS profile. For a ID Bjorken 
expansion and the SE profile Eq.(ll) gives a direct cor- 
relation between density p, temperature T and time r: 
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(12) 



with Tq,pq,tq the values at same initial time. By mean 
of Eq.(12) we can relate the AE/ At represented by the 
thick solid line in Fig. 1 to the delayed energy loss em- 
ployed by Pantuev [19]. In particular with our param- 
eters from Eq.(12) our thick solid line corresponds to a 
delay of about 1.8 fm close and even less extreme than the 
one in Ref. [19]. One can also notice that for some ob- 
servables already the less extreme GLV— a s (T) can give 
results similar to the low T energy loss, see Fig. 5 and 8. 

Hadronization - The final step of the model is the 
hadronization by independent fragmentation. The par- 
ton distribution after the jet quenching are employed 
to evaluate the hadron spectrum by independent jet 
fragmentation using the AKK fragmentation functions 
Dj(x,Q 2 ) which give the probability that a hadron H 
is formed from a parton of flavor /. The final hadronic 
spectrum is obtained from 



dN, 



H 



d 2 pj 



dx x 



E 



§±Df(x,Qi) (13) 



where x = pif jp T is the fraction of the / parton car- 
ried by the hadron H and Q — p T /2 is the energy 
scale. At RHIC there have been several evidence that 
while hadronization by independent fragmentation is able 
to describe proton-proton spectra at pt > 2 GeV, in 
Au+Au collision there are non-perturbative effects like 
quark coalescence that modify hadronization at least 
up to pt ~ 5 — 6 GeV [4, 5]. We do not include 
any hadronization by coalescence with the present work 
therefore all of the following results have to be considered 
reliable only for pt > 5 GeV. 



II. APPLICATION OF THE MODEL TO 
EVALUATE Raa{pt) 

The amount of quenching is quantified by comparison 
of the inclusive spectra d 2 N AA /dptdrj in ion-ion(AA) col- 
lision to a nucleon-nucleon(pp) reference cPcr NN jdp t dr\ 
via the Nuclear Modification Factor Raa{pi)- ■ 



RAAiPt) = 



d 2 N AA /dp t dr, 



T AA ■ d 2 a NN /d Pt dr 1 



(14) 



with Taa the nuclear overlap function which scales up 
single NN cross section to AA according to expected 
number of binary NN collisions without modification. 
Thus a Raa smaller(larger) than unity means suppres- 
sion(enhancement) due to medium effect. At RHIC this 
ratio at large pt > 6GeV has been measured to be nearly 
constant around avalue of 0.2 for the most central AuAu 



collisions, see Fig. 2 (upper). In our model the Raa(pt) 
can be calculated simply from the ratio of the spectra 
before and after quenching. 

We have applied our modelling of the jet quenching 
for Au+Au collisions at 200 AGeV. We use standard 
initial conditions that for the most central collision bin, 
b = 3/m, are given by a dN/dy = 1000 and a maximum 
temperature of the bulk To = 340 MeV at To = 0.6 fm, 
as usually done to describe the bulk in hydrodynamics 
and transport approaches [29, 30]. The results are shown 
in Fig. 2 and are performed for the two geometries de- 
scribed above: Wood-Saxon profile (WS) and Sharp El- 
liptic shape (SE). The GLV formula, Eq.(9), is used with 
a coupling a s — 0.27 but rescaled by a Z = 0.45 fac- 
tor [31] that accounts for higher order effect. Here it 
has been chosen to reproduce the data at pt = 6 GeV 
for the most central selection 0—5%. From the WS 
profile to the SE one has to decrease by about 15% the 
Z normalization factor due to lack of surface where the 
quenching is smaller. However both values are well in- 
side the uncertainty in the a s strength of the in-medium 
gluon radiation due to the uncertainty in the perturba- 
tive expansion and in the validity of the expansion itself. 
However as said in the introduction our purpose is not to 
constrain the total amount of quenching due to gluon ra- 
diation. Our methodology is to fit the R A a(pt) in order 
to fix the correct amount of total quenching with the aim 
of exploring the effect of different geometries and espe- 
cially different temperature (p or s) dependences of the 
energy loss on other observables. 

In Fig. 2 we can see that once the amount of quenching 
is fixed for the most central collisions both the depen- 
dence on pt and on centrality are correctly predicted with 
a GLV formula for quenching. Of course such a result 
has been obtained with other models of gluon radiation 
like BDMPS, ASW, AMY, DLGV [12-15]. Our purpose 
here was simply to show that our model is also able to 
reproduce the Raa(pt) which can be considered as the 
minimum requirement. However our result shows also 
that at the level of Raa(pt) both a realistic density pro- 
file like the WS or the Sharp Ellipse profile can describe 
the data reasonably well. We have seen that the main 
reason behind the similar final result between the two 
different geometries, WS and SE, relies on the compen- 
sation between two effects. In fact for SE one has a non 
realistic uniform density but on the other hand the mini- 
jet are also distributed uniformly. So the uniform den- 
sity profile from one hand overestimates the amount of 
quenching close to the surface but from the other under- 
estimates the one in the core of the fireball. In addition 
in such a modelling the fact that minijets are uniformly 
distributed overestimates the amount of minijets leaving 
the fireball nearly unquenched. Even if one could suspect 
that a balance among these effect should not apriori hold 
at all centralities our results show that the breaking of 
such a cancellation effects is small. In fact regard the 
Raa{pt) a simplified model with the SE cannot be dis- 
carded. The same conclusion can be drawn looking also 
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FIG. 2: Nuclear Modification Factor as a function of the 
transverse momentum px in Au + Au at 200 AGeV for differ- 
ent centralities. The circles are the experimental data taken 
from Ref. [11]. The calculation using the GLV energy loss 
are shown by dashed and solid line for a Sharp Ellipse pro- 
file and a Wood-Saxon profile respectively. The squares and 
open circles refer to GLV— a s (T) for a SE and WS profiles 
respectively. 



FIG. 3: Nuclear Modification Factor as a function of the num- 
ber of participant in Au + Au at 200 AGeV for two density 
profiles of the bulk matter: Glauber from Wood-Saxon nuclei 
profile (solid line) and Sharp Elliptic Shape (dashed line). 

pendence of the quenching even if one looks at the evo- 
lution with centrality. Therefore in agreement with Ref.s 
[32, 33], we find that Raa(pt) carries only a weak in- 
formation on the jet quenching process, apart of course 
the total amount of quenching, which in itself is of fun- 
damental importance and has lead to a first estimate of 
the average initial gluon density. 

In the next Section we will investigate both the elliptic 
flow at high p t and the flavor dependence of the quench- 
ing. In fact both observables are still hardly accounted 
for quantitatevely by the present models. 



III. ANGULAR AND FLAVOR DEPENDENCE 
OF THE QUENCHING 

The first analysis of jet suppression has shown that 
it is very difficult to have an agreement between mod- 
els and experiments for the dependence of the Raa on 
the azimuthal angle <j> respect to the reaction plane in 
non central collisions [33] . Such a dependence arise from 
the "almond" (elliptic) shape of the overlap region of 
two colliding nuclei. In particular for large pt > 6 GeV 
where hard processes dominate [11], partons penetrating 
the fireball in different directions lose different amount of 
energy according to their varying paths that on average 
are larger in the out of plane direction. A measure of this 
effect is provided by the second Fourier coefficient of the 
distribution, namely the elliptic flow: 



at the Raa at px > 6 GeV as a function of centrality ^ 

shown in Fig. 3. However, as expected, we see that sur- , . _ / # cos(20) [d 2 N / dp t d(f)} 

face effects leading to less suppression for a WS geometry v 2\Pt) — J 27T d(f> [d 2 N / dp t d(j)) 
are more important for a smaller N part . 

Furthermore if we consider a GLV formalism with Unexpectedly, measured V2{pt) happened to be con- 

a running coupling constant a s (T), GLV— a s (T), the siderably larger than what jet quenching models pre- 

Raa(pt) is well reproduced at same level of quality. dieted. It was noted in Ref.[34] that for very strong 

Hence looking at Raa{pt) one is not able to clearly dis- quenching only jets emitted from the surface of the "al- 

criminate neither the geometry nor the temperature de- mond shape" can be observed and the data for V2(pr) 



are very close to such a limiting case. However generally 
three other assumptions are made, namely: (i) quench- 
ing is proportional to matter density; (ii) colliding nuclei 
were approximated by homogeneous sharp-edge spheres, 
(iii) Only the net rate of energy loss is considered dis- 
carding the emission-absorption processes that generally 
can lead to an enhancement of the collective flows [13]. 
Studies by Drees, et al [35] relaxed the second assumption, 
with realistic nuclear shapes, which only made contradic- 
tion with data even stronger. A result that is confirmed 
also by the present work. 

We will mainly explore the effect of the first assump- 
tion by exploring the four different kinds of energy loss 
shown in Fig. 1. In Fig. 4, we show the resulting i>2(pt) 
for pions in Au + Au at ^s^n — 200 GeV and an impact 
pararameter b — 7.5 fm corresponding to a minimum 
bias condition for which the experimental value of ellip- 
tic flow is around 0.11 as shown in Fig. 10. We can see 
that even if the amount of total quenching has been fixed 
to the experimental value of Raa(pt) the amount of el- 
liptic flow is still strongly dependent on the temperature 
dependence of the AEi oss . Specifically from thin to thick 
solid line, we see that i>2 increases if the Ei oss is stronger 
as T — > T c (thick solid line of Fig.l) which means later 
in the evolution of the QGP, except for the surface. In 
Fig. 5 the same quantity is plotted for the case of a bulk 
density given by the Sharp Ellipse. Again we find that 
in agreement with other calculations [18, 19], the elliptic 
flow is significantly larger if most of the energy loss take 
place closer to T c (thick black line). If instead the en- 
ergy loss take place mainly at high temperature/density 
(thin black line) the V2(pt) is essentially vanishing. The 
dashed line is the result with the GLVc that represent a 
case in which the energy loss is essentially proportional 
to the density. The dot dashed line is the GLV-a s (T) 
that reweigths the temperature dependence through the 
coupling a s (T) giving an estimate of the effect of asymp- 
totic freedom. Such an effect is not negligible and for 
the SE geometry it seems to give already a V2{pt) very 
close to the much more extreme case rappresented by the 
thick solid line, i.e. quenching dominated by the close to 
T c region. 

The correlation between the temperature dependence 
of AE and the amount of elliptic flow developed is clear. 
In fact at variance with Raa(pt) the V2(pr) has a longer 
formation time because the minijets have to explore the 
shape of the fireball. For example if one consider the 
extreme case of a very strong quenching that takes place 
on distances much lower than the fireball size then there 
would be no elliptic flow (like in the Ei oss at high T case). 

It seems to be quite likely that experiments are telling 
us that quenching is likely to be not proportional to the 
matter entropy density or temperature, but a decreasing 
function of it. This is what is essentially discussed in 
Ref.s[18, 19] where however it was implicitely assumed 
that the amount of quenching of quarks and gluons are 
equal among them and to the hadronic one. Here we re- 
lax such assumptions showing that the temperature (or 
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FIG. 5: Elliptic Flow for the pions coming from quark and 
gluon fragmentation for the different T dependence of the 
energy loss as shown in Fig.l. The bulk density profile is a 
sharp uniform elliptic shape. 

entropy density) dependence of AEi oss modifies not only 
the V2(j>t) but also the relative amount of quenching of 
quarks and gluons. This does not want to be just a more 
detailed calculation but is indeed related to another puz- 
zle of the jet quenching phenomena observed more re- 
cently from the experimental study of the chemistry of 
the minijet suppression [36, 37]. We discuss it in the next 
subsection. 



A. Quark to Gluon Modification Factor 

The QCD due to its SU(3) Lie algebra gives a factor 
Ca = 9/4 larger for the energy loss of gluons respect to 
that of quarks. For this reason sometimes it is implicitely 
assumed that the ratio between the suppression of the 
gluons Raa{9) and the suppression of quarks Raa(i) is 
such that Raa(q) / Raa{q) = 9/4. From this expectation 
one would think that the (anti-) protons are more sup- 
pressed respect to pions because they come more from 
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gluon fragmentation than from quarks fragmentation re- 
spect to pions. The data at RHIC however has shown 
that also for such an observable there is no agreement 
with the data. In fact even outside the region where 
coalescence should be dominant [4, 5] the protons and 
the antiprotons appear to be less suppressed than the 
pions and p° [36, 37]. Again we can see that going be- 
yond the simple amount of quenching given by Raa(pt) 
neither the azimuthal dependence nor the flavor depen- 
dence of the quenching appear to be in agreement with 
the data. We call this open issues the "azimuthal" and 
the "flavor" puzzle respectively. We will show that even 
if Raa for central collisions is fixed to be ~ 0.2 also the 
Raa(q) I Raa[s) is affected by the temperature depen- 
dence of AEi oss 

In Fig. 6 we show the ratio of the Raa{<i) I 'Raa(9) for 
the WS geometry and the four different temperature de- 
pendences of the energy loss AEi oss as in Fig.l. We can 
see that the standard GLVc energy loss does not give the 
expected ratio 9/4 for Raa(q) / ' Raa{9) but a lower value 
around 1.8 which represents already a non negligible de- 
viation from 2.25. We can however see that if the energy 
loss would be strongly T dependent and dominated by 
the T > T c region Raa(q) / Raa{9) can increase up to 
about 2.3 while oppositely if it is dominated by the high 
temperature region (thin solid line) the Raa{9) / 'Raa(9) 
can become as small as 1.5. On the other hand our study 
of the elliptic flow, as well as previous studies, show that 
an energy loss that increase with T would generate a tiny 
v%{pt) very far from the observations. Such an effect is 
totally discarded in ref.s[18, 19] that neglect the different 
quenching of quarks and gluons. In other words while the 
indication of a AE(T) increasing as T — > T c is confirmed 
also by our model in order to reproduce experimental 
data for vi , we have spot that this would lead to a larger 
Raa(q) I ' Raa{9) ratio increasing the disagreement with 
the Raa(pt) observed for the various hadrons like p,p, n 
that show the Raa of baryons less suppressed respect to 
the pionic one. 

In this respect we spot also another potential prob- 
lem that can arise for simplified fireball modelling when 
extreme AE(T) are considered and more exclusive ob- 
servables are investigated. More explicitly we refer to the 
Sharp Ellipse (SE) fireball that neglects the surface of the 
fireball. In Fig. 7 we show that the results have an oppo- 
site pattern for the SE geometry respect to the realistic 
WS one. The standard GLV energy loss gives approx- 
imately the expected ratio of 9/4 for Raa{q) I Raa(9)', 
This probably led to associate the 9/4 factor on Ei oss 
to that on Raa(q) / Raa(9) ■ However this ratio is more 
sensitive to the temperature dependence of the energy 
loss, in a way that is exactly opposed to that we have 
seen for the WS geometry. For this kind of geometry 
(SE) if the energy loss is stronger at high temperature, 
the Raa(q) / Raa(9) can result to be almost 5 (thin solid 
line) on the other hand, if the enery loss is stronger closer 
to T c the Raa would be just sligthly above 1 (thick solid 
line). This means that the ratio tends to decrease from 
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dence of the energy loss as shown in Fig.l. The short-dashed 
line corresponds to the 9/4 value. The density profile of the 
bulk is that of a sharp elliptic shape. 



5 to 1 if the energy loss is dominated by the high tem- 
peratures (T ~ 2T C ) or by the low ones (T ~ T c ). 

It is instructive to explain the origin of such a strong 
effect of the density profile, but before we note that such 
an un-expected strong effect shows up only when extreme 
AE{T) are considered. In fact the standard GLVc Ei oss 
is modifed by about a 15% moving from a WS-Glauber 
to a simple SE profile. 

In order to understand the mechanism behind the de- 
termination of the ratio Raa{q) / Raa(.9) we discuss an 
oversimplified case in which all quarks lose the same 
amount of energy and all gluons lose their energy ac- 
cording to AE g — %/AAE q . For such a simple case the 
spectra after quenching are shifted by a quantity equal 
to the lost energy. Quarks that finally emerge with an 
energy Ef = px are those which before quenching had an 
energy E t — p T + AE. So Raa for quarks is equal to the 
ratio between the parton distributions in momenta space 
without quenching J(pt) and the quenched one given by 
f(pr + r] AE) where f{pr) = dN/ d 2 pr dy is given by Eq. 
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(1), rj = 1 for quarks and rj = 9/4 for gluons. Therefore 
the Raa is 



Raa(pi 



f(PT + V AE) 

KPT) 



(16) 



and the ratio between quark and gluon nuclear modifica- 
tion factors is 



R A A{q) _ f q (p T + AE) 



f g (PT) 



(17) 

RAA(g) f q {p T ) fg (j>t + (9/4) AE) v ' 

of course there is no reason why this ratio must be 
9/4 and we can also see that even without the 9/4 factor 
there can be a Raa(q) /Raa (g) that is not one. In Fig. 8 
is shown the dependence of the ratio on AE for partons 
with p T = 10 GeV and we can see that Raa(q) I Raa{,9) 
quickly increase with AE even if the ratio between the 
quark to gluon energy loss is 9/4. The ratio is about 4 for 
a AE ~ 2 GeV/fm. We will see that this observation is 
fundamental to understand the relation between AEi oss 
and the ratio Raa{q) / Raa{q) ■ In this very simplified 
model the Raa(q) / RAA{g) would be very large for an 
average energy loss < AE q >= AGeV typical of central 
collisions where Raa(pt) ~ 0.2. From Fig. 8 we can see 
that this would give a Raa{q) / Raa(q) ~ 8. 

One can move toward a minimal realistic model dis- 
tinguishing among partons two classes of particles: those 
that undergo a large quenching and those that lose no 
energy or better a small amount of energy, usually asso- 
ciated with the minijets generated at the surface. Let's 
consider a case in which 50% partons lose an energy 
AE q = lGeV = 4/9 AE g while the other 50% lose the 
energy AE q indicated in Fig. 8. In such a case the evo- 
lution of Raa(<i) I Raa(9) is very different respect to the 
first case as shown by the much milder increase of the 
dashed line that reaches a maximum of about 2.6 and 
even decrease for a AE q > 2 GeV. 

This clearly indicates that because of the rapidly 
falling distribution with px, the determination of the ra- 
tio is dominated by those partons that suffer less energy 



loss . If one for example suppose that 50% of the particles 
do not lose energy then whatever is the energy loss of the 
other partons the Raa(q) / RAA(g) will be equal to one. 
This is shown by the thin solid line in Fig. 8. Indeed al- 
ready if only a 10% of partons do not lose energy the ratio 
stays below 9/4 whatever is the Ei oss of the other 90% of 
particles and even decrease if those particle lose a large 
amount of energy bringing the value again closer to one. 
This is shown by the dot dashed line in Fig. 8 where one 
can realize the huge impact of the particles that do not 
lose energy comparing thick solid and dot-dashed lines 
which differ only by the fact that a 10% particles do not 
lose energy. Finally with the dot line in Fig. 8, we show a 
system in which most of the particles undergo a quench- 
ing with AE q = 2.5 GeV = 4/9AE g . This is closer to 
a case in which most of the particles undergo a simi- 
lar strong quenching. In summary we can say that the 
Raa(q) I Raa{s) is not really directly linked to the rela- 
tive amount of quark and gluon energy loss because it is 
largely affected by the way the energy loss is distributed 
among partons. In particular once there are minijets that 
do not suffer energy loss the ratio Raa(q) / RAA(g) gets 
closer to one, because it is more affected by these mini- 
jets. Hence also a careful treatment of the corona effect 
would much likely give a significant contribution to the 
determination of the Raa{q) / 'Raa{.9) ■ 

On the base of the above discussion it is possible to 
understand the dependencce of R AA {q) I RAA(g) on the 
temperature dependence AE(T), seen in Fig. 6 and 7 and 
its opposite behavior between the WS and SE geometries 
for the density profile. In the case of energy loss at low 
temperature with a SE profile the quenching happens 
at the end of the lifetime of the fireball because there 
is a direct relation between time and temperature, see 
Eq.(12). Therefore in the case of the energy loss increas- 
ing as T —¥ T c many particles escape without losing en- 
ergy for SE profile. Only the particles in the inner part 
of the fireball are quenched. This means that we are 
closer to the schematic case described by the thin solid 
line for large AE in Fig. 8, most particles do not lose 
energy the rest lose a large amount of energy, and in fact 
Raa(q) I Raa{<3) is close to one, see Fig. 7. 

Instead if the quenching is larger at high temperature 
all particles lose energy early, and except for those very 
close to the surface most of particles AE > 2.5 GeV. In 
this case there are essentially no particles that do not 
lose energy. So we are closer to the case described by the 
dotted line for which Raa(q) / RAA(g) is about 5 to be 
compared to thin solid line of Fig. 7. 

In the case of WS geometry in general there is no di- 
rect relation between time and temperature because on 
the surface one will have already at early times a low 
temperature. Therefore in this case an energy loss dom- 
inated by high temperatures means that the quenching 
is large only on the inner part of the fireball while par- 
ticles in the surface lose a small amount of energy. We 
are close to the case described by the dot dashed line in 
Fig. 8 and the R AA (q) / R AA (g) is in fact about 1.5. It is 
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difficult to reach one because due to the density profile 
anyway all the particles lose at least some finite energy. 
Instead if the quenching takes place mainly at low tem- 
perature (T ~ T c ), this means that energy loss is strong 
in a layer on the surface of the fireball, and because all 
particles must go through this layer at some time, all of 
them lose a large amount of energy. We understand that 
this is very different, essentially the opposite, respect to 
SE case. In the latter case a AE(T) with a maximum 
at T ~ T c means that most of the particles escape due 
to the direct time-temperature relation that is not de- 
pendent on space. For WS-Glauber profile such a case 
means just the opposite, i.e. all the particles lose a similar 
amount of energy. This is essentially the reason under- 
lying the opposite trend seen in Fig. 6 and 7. We finally 
notice that the importance of such details emerge only if 
extreme energy losses T-dependence are considered. 

The conclusion of this study is twofold: i) If one tries 
to reproduce the large value of elliptic flow using a type 
of energy loss that increase with decreasing temperature 
there is a simultaneously increase of the Raa(q) / 'RAA(g) 
enhancing the discrepancy respect to the observed "fla- 
vor" dependence of the suppression that seems to prefer 
a ratio close or even smaller than one [36] . 

ii) When peculiar energy dependences are considered 
the specific density-temperature profile can become very 
important for a quantitative evaluation of the observ- 
ables. 

These results bring us to make two important obser- 
vations: 1) The SE profile is able to describe the ob- 
served Raa(pt), but it is inadeguate to reproduce the ra- 
tio Raa(q) I Raa{s) j and furthermore, it cannot be used 
even for a qualitative analisis of this ratio because it gives 
opposite results to the more realistic WS profile. 2) In 
Fig.s 4 and 5 we have shown that an enhancement of 
the energy loss near T c increases the elliptic flow toward 
an agreement with the experimental data. On the other 
hand such a AE(T) is associated to an enhancement of 
Raa(q) / RAA(g) in apparent disagreement with the data 
[371. 



IV. CORRELATION BETWEEN Raa (q)/ Raa {g) 
AND ELLIPTIC FLOW 



To solve what we have called the " flavor puzzle" inelas- 
tic collisions that cause a change of the flavor has been 
invoked [20-23] . Such a process would at the end produce 
a net conversion of quarks into gluons. Hence a decrease 
of gluon suppression respect to the direct suppression 
and an increase of the quark one. [20]. In Ref. [20] it 
has been calculated the conversion rate of a quark jet 
to a gluon jet and vice versa due to two-body scatter- 
ings. An enhancement factor K c = 4 — 6 that accounts 
for non-perturbative effect is needed to produce a nearly 
equal suppression of quarks and gluons. This is not an 
unreasonable enhancement factor considering that at our 
energy a K ~ 4 is necessary also to have the right mini- 
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FIG. 9: Ratio of quark to gluon in Au + Au collisions (b=7 
fm) for the energy loss derived in Ref. [20] with different factor 
upscaling the pQCD conversion rate. 



jet initial distributions in pp if one starts from a simple 
second order pQCD calculations. We have included the 
flavor conversion mechanism in our code in a fashion sim- 
ilar to Ref. [20]. In order to check our code and to have 
a direct link to the previous calculations of Fries, Ko and 
Liu we have used the same density profile (i.e. a SE pro- 
file), same conversion rate and AEi oss derived by them 
at leading order in the pQCD expansion. 

In Fig. 9 we show Raa{<1) I Raa{9) for different values 
of K c . Notice that with the AEi oss of Ref. [20] the ratio 
Raa(q) I ' Raa{9) without conversion is again different and 
moreover quite larger than 9/4. However a K c ~ 4 — 6 is 
able to reduce that ratio by about a factor of three mak- 
ing it close to unity. After this test we have fixed K c = 6 
and performed the calculation with the four different en- 
ergy loss as shown in Fig.l. In Fig. 10 we show directly 
the (Raa{q) / Raa(q) , V2) plot. This plot manifest a clear 
but non-trivial correlation between Raa(q) / Raa(9) and 
the V2(pt)- The upper symbols are the results without 
the jet flavor conversion and we can see that such a cor- 
relation drives the data far from the experimental ob- 
served values of a t>2 ~ 0.1 and an Raa {<i)IRaa (5) _ 1 
(to account for the Raa{p + P) > Raa(^ + + tt~) with 
AKK fragmentation function). The lower symbols corre- 
sponds to the results including the rate of inelastic col- 
lisions and we can see that this process allows to get 
closer to the experimental region because they change the 
Raa{q) I Raa(9) without modifying the elliptic flow of pi- 
ons. Fig. 10 essentially demonstrates that if one tries to 
reproduce simultaneously the Raa{pt) and ^(pt) and 
Raa{q) I RAA(g) a AE(T) increasing as T — > T c is needed 
but also that in such a case a flavor conversion process 
becomes even more necessary. 

However at this point a quantitative study should be 
performed employing a more accurate dynamics for the 
bulk evolution given by hydrodynamics or partonic trans- 
port theory that include three-body radiative processes 
[41, 42]. In fact, especially if more peculiar energy loss 
dependence has to be investigated, it is important to have 
a quite realistic density and energy density profile as we 
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have discussed above. Furthermore fluctuations in the 
energy loss [10, 38], the gain-loss processes [13] and the 
elatic energy loss [39, 40] should be included because they 
can give correction to both Raa(pt) and V2(pr)- 

V. IMPACT OF THE EQUATION OF STATE 

As a last point, we explored the impact of the equation 
of state (EoS) on the observables and especially the cor- 
relation between Raa (<7) / Raa (<?) an d V2{pt)- This has 
to be taken as a first explorative study. 

We notice that even if the free gas approximation is 
a quite reasonable approximation to describe the rela- 
tion between density and temperature for most of the 
evolution of the expanding QGP it is not true close to 
T c where the relation between temperature and density 
is strongly modified by the cross-over phase transition. 
One generally thinks that for the description of the high 
Pt particles this can be discarded in first approximation. 
However once it is opened the way to the possibility of a 
quenching that is not proportional to the density but is 
stronger close to the phase transition we will show that 
one should more carefully look into the problem. In fact if 
quenching would be dominated by the T ~ T c region the 
time spent in this region by the fireball can be strongly 
modified if a realistic EoS is considered. In fact for a 
simple 1+1D expansion the relation between tempera- 
ture and density is modified from Eq.(12) to: 

^=(^y (T, =P)-> (18 , 

To \PoJ W/ 

where P{T) is a temperature dependent coefficient that 
can be obtained from a fit to lattice QCD data [43]. We 



have found: 

f3(T) = \-a(^y , T>T C (19) 

with a = 0.15 and n = 1.89 and of course for T >> T c 
one gets (3 ~ 1/3. The parameters has been calculated 
from a fit to the lattice QCD data of Ref.[43] on the 
energy density and pressure. We see that close to T c the 
/3 coefficient is quite small which means that even if the 
density goes down as r _1 , the temperature stays nearly 
constant as one can expect in a first order or strong cross- 
over phase transition. This means that the system spends 
more time close to T c respect to the simple picture given 
by Eq.(12) usually assumed. 

In order to estimate the impact of this correction we 
have performed a simulation for the AEi oss (T) behavior 
represented by the thick solid line in Fig.l which is similar 
to the delayed energy loss proposed by Pantuev [19] as a 
solution for the observed large elliptic flow. We consider 
only this case because it is of course the one that is much 
more affected by the modification implied by Eq.(18). 
Instead for the opposite case of an Ei oss dominated by 
the high T the effect is vanishing. We have again taken 
care to regulate the total energy loss in order to have the 
correct amount of Raa(pt)- 

The results are given by open symbols in Fig. 10. The 
open square is the result without the inclusion of the fla- 
vor conversion mechanism and shows a further increase 
of both V2(pt) and Raa(q) I 'Raa{9) , in line with the 
previously seen correlation. The enhancement of the 
V2{pt) in itself is due to the fact that with a EoS the 
system spend more time close to T c , therefore in or- 
der to have the same amount of suppression most of 
the energy loss occurs later. However more interest- 
ingly the inclusion of the q O g conversion results to 
be more efficient that in the previous case generating 
a strong decrease of Raa{<i) / 'RAA(g) while keeping the 
same V2{pt)- The stronger effect of flavor conversion 
is again due to the longer lifetime. The final result is 
a combination that move the final value of both V2(pt) 
and Raa(<1) / Raa(9) much close to the experimental ones 
shown by the shaded area in Fig. 10. We also notice that 
a Raa (9) I Raa (g) < 1 is obtained while generally it is 
believed that flavor conversion by inelastic collisions can 
give at most Raa(q) / RAA(g) = 1 [36]. 

This last result about the impact of a realistic EoS de- 
serves a more careful treatment again employing a more 
realistic description of the bulk. Indeed a full description 
of the dynamics related to the cross-over region should in- 
clude also the gradual change from a quark-gluon plasma 
to an hadronic gas. However this would be out of our aim 
for the present work. 

VI. CONCLUSIONS 

In the present work we have considered different tem- 
perature dependences of the energy loss tuning always 
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the parameters to reproduce the experimentally Raa(pt) 
suprcssions of pions. We have shown that even if the 
Raa(pt) is fixed, different AE(T) generate very differ- 
ent value for both the V2(pt) and the Raa(q) I 'Raa(9)- 
Indeed both quantities are quite puzzling because stan- 
dard jet quenching models are not able to reproduce none 
of them. We refered to this as the "azimuthal" and 
"flavour" puzzle. We have found however that there is 
a correlation between these two observables that is de- 
tcmined by the temperature (or density) dependence of 
the quenching. In agreement with Ref.[18, 19] we have 
found that if the quenching is dominant closer to T c the 
V2(pt) is enhanced getting closer to the data. However 
while in Ref.[18, 19] it is discarded a separate treat- 
ment of quark and gluons, their Ei oss leads to a quite 
large ratio RAA{q)/ Raa{9) which would result incom- 
patible with the observed systematics of Raa(p + p) > 
Raa{k + + tt~) ~ Raa(p°)- It appears that while the 
V2(pr) would suggest a AE(T) that increases towards T c 
the Raa {q) I Raa (g) would become too large even if we 
do not yet have a direct measurement of Raa{q) / RAA(g) 
and we impinge on the uncertainties coming from the 
fragmentation functions. 

In this context we have also spot an unexpected strong 
dependence on the density profile of the fireball that 
emerge for extreme choices of AE(T), while the depen- 
dence is milder but not negligible for more standard en- 
ergy loss like the GLV one. This puts a warning for 
further studies: once going beyond the simple Raa(pt) 
factor one has to rely on a realistic dynamical evolution 
of the bulk matter like those supplied by hydrodynamics 
and/or parton cascades [41]. 

It seems that the only way to get closer to the ob- 
served values for R A a(pt), v 2 {pt) and Raa{i) / Raa{9) 



solving both the "azimuthal" and "flavour" puzzles, is 
to have both AE(T) increasing close to T c and a flavour 
conversion mechanism. Finally we point out that while 
generally the free gas approximation is quite reasonable 
to describe the expansion of the QGP fireball, this is no 
longer true if the energy loss occurs dominantly close to 
T c . In such a case it has a significant impact to take into 
account the strong deviation from the free gas approxi- 
mation occuring in the cross-over region and leading to 
a slowing down of the cooling close to T c . This increase 
further the time spent around T c which enhance both 
V2(pt) and the efficiency of the q <R- g in medium conver- 
sion moving the values of (Raa(q) / RAA{g) , V2(pt)) close 
to the observed ones. 

It is of course important to study how the longer life- 
time and higher temperatures reached at LHC energies 
affects the observed correlations. Furthermore in Ref.[24] 
it has been pointed out that a better probe of the fla- 
vor conversion mechanism is supplied by an high px 
strangeness enhancement. It remains to be studied if 
such an enhancement is affected by the T dependence of 
the energy loss. 

In Summary we have pointed out the impact of pe- 
culiar T dependences of the Ei oss on both v 2 (pt) and 
Raa (9) I Raa (g) and their correlation. Moreover we have 
spot the relevance that the EoS may have in case of 
Ei oss dominated by the T ~ T c region. In any case our 
study, although already revealing several interesting indi- 
cations, is mainly explorative. A more quantitative anal- 
ysis should be performed with more sophisticated models 
that include the energy loss fluctuations, realistic gain 
and loss processes, elastic energy loss and a more accu- 
rate description of the bulk. 
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